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Abstract 

We prove that arbitrary (non-polynomial) scalar evolution equations of order m > 7, that 
are integrable in the sense of admitting the canonical conserved densities p^^\ P^^\ and p'^^^ 
introduced in [A.V. Mikhailov, A.B. Shabat and V.V Sokolov "The symmetry approach to 
the classification of integrable equations" in ' What is Integrability? edited by V.E. Zakharov 
(Springer- Verlag, Berlin 1991)], arc polynomial in the derivatives u^-i for i ~ 0,1,2. We 
also introduce a grading in the algebra of polynomials in Uk with fc > m — 2 over the ring of 
functions in x,t,u,.. . , Um-s and show that integrable equations are scale homogeneous with 
respect to this grading. 
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1. INTRODUCTION 



The classification problem for scalar integrable equations in one space dimension is solved in the 
work of Wang and Sanders [1] for the polynomial and scale invariant case where it is shown that 
integrable equations of order greater than or equal to seven are symmetries of third and fifth 
order equations. In subsequent papers these results are extended to the cases where negative 
powers are involved [2] but their methods were not applied to equations without polynomiality or 
scaling properties. 

The aim of the present work is the study of the classification problem for arbitrary evolution 
equations. In our work, we use the "formal symmetry" method, where the existence of certain 
conserved densities is a necessary condition for integrability. Our main result is that evolution 
equations that are integrable in the sense above are polynomial in top three derivatives and 
have a certain scale homogeneity property. 

The "formal symmetry" method introduced by Mikhailov, Shabat and Sokolov [3] have 
been used to obtain a preliminary classification of essentially non-linear third order equations 
and quasi-linear fifth order equations. The classification of essentially non-linear equations is 
considered later in the work of Svinolupov [4], where quasi-linear integrable equations that 
are not linearizable are found to be related to the Korteweg-deVries and Krichever-Novikov 
equations through differential substitutions. 

The first result towards a classification for arbitrary m'th order evolution equations is 
obtained in [5] where it is shown that scalar evolution equations ut = F[u], of order m = 2k + l 
with m> 7, admitting a nontrivial conserved density p = Pu^ + Qun + R of order n = m+1, are 
quasi-linear. The method of [5] is not applicable to third order equations, because for m = 3, 
the canonical conserved density p^^^ is not of the generic form on which the quasi-linearity result 
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is based on. For m = 5, although the generic form of p^^) is vahd, k = 2 occurs as an exception 
and we cannot exclude the existence of fifth order non-quasi-linear integrable equations, at least 
with the present method. For m > 7 the structure of integrable equations seem to be much 
simpler and one may hope to obtain a complete classification as in the polynomial case. 

In the present paper we continue with the classification problem using formal symmetries 
and we prove that evolution equations of order m = 2k + 1 > 7 admitting the canonical 
densities pW, i = 1,2,3, as given in Appendix A, are polynomial in the derivatives Ufni ^^m— i 
and Um-2- The final result presented in Corollary 4.6.2 gives the explicit form of a candidate 
for an integrable evolution equation as a polynomial in Um, Um-i and Um-2, with coefficients 
as yet undetermined functions of lower order derivatives. This result is definitely the best one 
can obtain by the use of the canonical densities p^'^\ for i up to 3, as discussed Remark 4.7 and 
in the conclusion. Any further progress towards polynomiality would require the computation 
of new canonical densities. These computations for general m are extremely tedious and they 
are deferred to future work. There is also an alternative and more promising direction towards 
the classification: Our candidates for integrable equations have a certain scale homogeneity 
property with respect to a grading in the algebra of polynomials in Uk with k > m — 2 over the 
ring of functions in x,t,u, . . . ,Ujn-3- This grading called the "level grading" and introduced 
in [8] proved to be an efficient tool for practical computations and it is expected to allow the 
treatment of the polynomiality in lower orders in a unified manner. However, as both directions 
of approach to the problem require completely different techniques, the present paper is limited 
to the information that can be extracted from the existence of the conserved densities p^*) , for 
i < 3. 

The notation and terminology are reviewed in Section 2 and classification results are given 
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in Sections 3 and 4, where we show that for sufficiently large m, integrable equations of order 
m are polynomial in the top 3 derivatives, Um, Um-i and Um-2- These results are obtained 
from the requirement that the canonical densities p^'^\ i = 1,2,3, be conserved quantities. The 
expression of these canonical densities are given in Appendix A. In Appendix B, we present 
general formulas for the k'th order derivatives of differential functions for large k. Due to 
the restrictions on the validity of these derivatives, the general formulas obtained in Section 
4 are valid for m > 19. For m < 19, we have done explicit computations with the symbolic 
programming language REDUCE and obtained the same polynomiality results, as summarized 
in Section 5. 



2. NOTATION AND TERMINOLOGY 

Let u = u{x, t). A function ip of x, t, u and the derivatives of u up to a fixed but finite order will 
be called a "differential function" [6] and denoted by (p[u\. We shall assume that <^ has partial 
derivatives of all orders. We shall denote indices by subscripts or superscripts in parentheses 
such as in a(j) or p^^ and reserve subscripts without parentheses for partial derivatives. For 

example if u = u(x,t), then 



du 



u 



while for ip = ip{x, t,u,ui, . . . , Un), 

dip dip dip 

For (p as above, the total derivative with respect to x is denoted by Dip and it is given by 



Dip = ^ipiUi+i + ipx. (2.1) 

i=0 



Higher order derivatives can be computed by applying the binomial formula 

'k-l 



1=0 



3=0 
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+ D^-\^. (2.2) 



The total derivative with respect to time denoted by Dt is given by 

n 

Dt<f = J2'PiD'F + <ff (2.3) 

i=0 

We recall that a differential function p is called a conserved density, provided that there is 
a differential function a such that Dtp = Da. If p is polynomial in certain higher derivatives, 
in order to check the conserved density condition, one can proceed with integration by parts 
and require the vanishing of the terms that are nonlinear in these highest derivatives. In our 
derivations we shall use only the vanishing of the coefficients of top two nonlinear ities. As 
shown in Proposition 3.1, these terms come from top 4 derivatives in the expansion of Dtp. 
The general expression for D^tp given by {A.6d) is valid for A; > 7 and as we use the derivative 
£)k-2p ^Yie general formulas, we need A; > 9 hence m > 19 for the validity of the general 
expressions. 

We shall denote generic functions (p that depend on at most by 0{un) or by \(p\ = n. 
That is 

(f = 0{un) or l^^l = n if and only if — = for A; > 1. 

OUn+k 

li if = 0{un), then Dip is linear in Un+i and D'^ip is polynomial in Un+i for z > 1. In order to 
distinguish polynomial functions we use the notation ip = P{un), i.e., 

p = P{un) if and only if ip = 0{un) and — — r- = for some k. 

This distinction have been used in the expression of derivatives given in Appendix B. 

Note that the total derivative with respect to x increases the order by one, thus if \p\ = n 
then l-D^f/jl = n + k. Furthermore, when ut = F, with \F\ = m, Dt increases the order by m. 



Equalities up to total derivatives with respect to x will be denoted by =, i.e., 

ip = tp if and only if tp = ip + Drj 

We shall repeatedly use integration by parts of the following type of expressions. 
Let pi < p2 < ■ ■ ■ < pi < s — 1 and \(p\ = k < pi. Then 

These integrations by parts are repeated until one encounters a monomial which is nonlinear 
in the highest derivative. 

The order of a differential monomial is not invariant under integration by parts, but one 
can easily compute the order of the nonlinear term that will result after integrations by parts. 

3. GENERAL RESULTS ON CLASSIFICATION 

We start by a brief description of the formal symmetry method [3] . If i? is a recursion operator 
for the evolution equation ut = F[u], then any fractional power of R is also a recursion operator. 
Thus starting from a recursion operator of order 1, expanded in a formal series in inverse powers 
of D, one can compute the operators of orders k, for k = 2,3, . . .. Each of these operators 
satisfy the operator equation 

R'^ + [R\F4 = 0, 

where is the Frechet derivative of F, defined by = YliLo FiD'^. It is known that for any 
two formal series A and B, the coefficient of is a total derivative [7], hence the coefficient 
of in each R^, denoted by p^^\ is a conserved quantity. In addition to these, for any m, it 
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is known that 

are conserved densities. We recall that as stated in Section 2, subscripts denote differentiation 
with respect to the derivatives of u, that is = and pm,n denotes ^ , etc. 

The explicit expressions of p^^^ and p^^^ for m > 5 and of p^^^ for m > 7 obtained in [5] are 
given in Appendix A. 

In order to compute Dtp up to total derivatives, we use (2.3) and integrate by parts 
until we encounter a term which is nonlinear in the highest derivative. The derivation of 
the polynomiality result is based on the coefficients of top two nonlinear terms, given by the 
equations (3.2) and (3.3), in Proposition 3.3. To arrive to these expressions we first show in 
Proposition 3.1 that if F has order m = 2k + 1 and p has order n = m + l, then, for m > 19, 
the top two nonlinear terms are u'^j^_^i_^_i and u^j^_^i- 

Proposition 3.1 Let p = p{x,t,u, . . . ,Un) be a differential function of order n and ut = 
F(x,t,u, . . . ,Um) be an evolution equation of order m, where m = 2k + 1, n = m + I and 
A; + Z — 1 > 0. Then, up to total derivatives Dtp is 

{-l)'+'Dtp - (D'^+'pn - D^pn-l) D^+'F - (D^Pn-2 - D'^-'pn-s) D^+'-^F 

(3.1) 

where ip = 0{u-ik+i-i)- 

Proof. In Dtp = YA=oPi^^^ + Pti highest order derivative comes from p^iD^F, where pn 
and D'^F are of orders 2k + I + 1 and Ak + I + 2 respectively. If we integrate by parts k + 1 
times we obtain 

PnD^'F ^ {-l)''+^D''+^pn D^+^F, 
6 



where D^^^pn and D'^'^^F are now respectively of orders 3k + l + 2 and 3A; + Z + 1. One more 
integration by parts gives a term nonlinear in usk+i+i- Similarly one can see that in 
Pn-i and D"^~^F are of orders 2k + l + l and Ak + l + 1. This time, integrating by parts k times, 
we have 

where D^pn-\ and D^^^F are both of orders 3fc + Z + 1. Thus the highest order nonlinear 
term in u:ik+l+i comes from top two derivatives in pnF>"' and pn-\D^~^F. By similar counting 
arguments, one can easily see that top two nonlinear terms are obtained from top four derivatives 
and the remaining terms are of order 3fc + Z — 1. □ 

Remark 3.2 As the general expressions for the derivatives given in (A.6a-d) are valid for large 
k, there are restrictions on the validity of the formula (3.1). Since the top four terms of D^^^F 
and D^^^pn are needed in (3.1), from (A.6d) it follows that k+ 1 and k + l should be both larger 
than or equal to 7. On the other hand, at most top two terms of the expressions in the second 
bracket in (3.1) contribute to the top nonlinearities and it turns out that the restrictions coming 
from (A.6a,b) are always satisfied and the crucial restriction is + 1 > 7 and A; + Z > 7. Thus 
for Z = 1, 0, —1 and —2, we need respectively k > 6(m > 13), k > 7{m > 15), k > 8(m > 17) 
and k > 9(m > 19). 

We shall now give the explicit expressions of the coefficients of top two nonlinear terms for 
m > 19 and k + l>7. 

Proposition 3.3 Let p = p{x,t,u, . . . ,Un) be a differential function of order n and ut = 
F{x, t,u, . . . , Ujn) he an evolution equation of order m, where m = 2k + l,n = m + I. Then for 
k > 9 and k + I > 7, the coefficients of the top two nonlinear terms u1|^_^_l_^_^ and u^|^_^_l in the 
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expression of Dtp up to total derivatives is 



{-if+^Dtp = Ai ulk+i+i + Ao ulk+i + V 



where ip = 0(u3k+i-i) and Ai and Aq given below. 



Al = {k + \)FmDpn,n - {k + I + \)DFmPn,n - Fm-lPn,r, 



Ai = Pn,n D^Fm Ijo C^k^ + 6A;^^ + 6kl'^ + 21^ + 3A;^ + + 6kl + k + l) 



+ pn,n D^Fm-1 \ {k^ + 2kl + 2k + 21 + f + I) 



+ pn,n DFm-2 \\ (3 + 2A; + 20 



"I" Pn.n Fjfi—3 



+ £)p„,„ D^Fm \l {-2k^ - Ak^l - 2kf + k^ + 1^ + 2kl + k + I) 



+ Dpn,n DFm-1 \ {I + I - 2k'^ - 2kl) 



+ Dpn,n Fm-2 \\ (1 - 2k) 



+ D''pn,nDFm \{2k'^ +2kH-k'' -k) 



+ D^p^^ri Fm-1 [I k^ 

+ D^pn,n F^ (-2fc3 - - k) 

+ Dpn,n-1 DFm [| (-1 + 2k + 21) 

+ £>Vn,n-l Fm [\ ("1 " 2fe) 

+ Pn,n-2 DFm [^k + 21 - \] 



+ Pn,n-2 Fm-\ [2] 



+ Dpn,n-2Fm [-2A; - 1] 

+ Pn-l,n-l DFm [\ {I - 2k - 21) 

+ Pn-l,n-l Fm-1 [—1] 

+ Dpn-l,n-l Fm \\ (1 + 2A;) 



(3.3) 



Proof. The proof is a straightforward computation of the integrations indicated in Proposition 
3.1. Writing the first four terms in Dtp and keeping only the terms which contribute to the 
nonlinearities , we get 

{ — l)^^^DtP = Pn,n FmUsk+l+1 U3k+l+2 

+ Pn,n [Fm-1 + {k + l)DFm] U3k+l «3fc+Z+2 



+ 
+ 



+ 



+ Pn,n FjYi—2 

+ ik + l)DFm-i + {''t)D^F, 



U3k+l-l U3k+l+2 



+ pn,n [Fm-3 + {k + l)DFm-2 + {^tW'Fm-l + CtW^F, 
+ (k + l) Dpn,n Fm U^k+l+l '"3fc+Z+l 



+ {k + l) Dpn,n [Fm-1 + {k + l)DF^] U3k+l U3k+l+l 



'"3fc+i-2 U3k+l+2 



+ ik + l)Dpn,n Fm-2 + {k + l)DFm-l + {^f)D^F, 



UZk+l-1 U3k+l+l 



Pn,n-2 + Dpn,n-1 'Pn,n — Pn-l,n-l FmU^k+l+l U^k+l 

Pn,n-2 + Dpn^n-1 + C2^)F>'^Pn,n — Pn-l,n-l 

X [Fm-1 + {k + l)DFm] U3k+l U^k+l 
Pn,n-3 + {k + l)Dpn,n-2 + kD^ Pn,n-1 + Ct^)F>^Pn,n - Pn-l,n-2 - kDpn-l,n-l 

xFm U3k+l+l U3k+l-l 
Pn-2,n Fm U3k+l ^3fe+/+l 



- Pn-2,n [Fm-1 + {k + I - l)DFm] U^k+l-l U^k+l+l 



- [Pn-2,n-l + kDpn-2 ,n Pn—3,n\ -^m U3k+l U3k+l- (3.4) 

After integrations by parts we get (3.2) as the coefHcient of the first nonlinear term 

and (3.3) as the coefficient of the second nonhnear term □ 

Actually the expression of Ai is valid for m > 13 and A; + / > 7. It is clear that if p is a 
conserved density for the evolution equation ut = F, then necessarily one should have 

Ai = 0, Ao = 0. 

These two equations will be used repeatedly in order to derive a number of necessary condition 
for integrability. 

Prom equation (3.2) we can easily get a number of results pertaining the form of the 
conserved densities. In particular we can see that higher order conserved densities should 
be quadratic in the highest derivative and top coefficients of the conserved densities at every 
order are proportional to each other [5] . 

Corollary 3.4 Let p = p{x,t,u, . . . ,Un) be a differential function of order n and ut = 
F(x, t,u, . . . , Urn), be an evolution equation of order m, m> 7 and n > m. Then 

Pn,n,n = (3.5) 

Proof. It can be seen that (3.2) uses only top two terms and for i > it is valid for + 1 > 3. 
Writing it in the form 

ik + + l + = (3-6) 

Pn,n 

we can see that for n > m the highest order term is Dpn^n and it follows that Pn,n,n = 0. 1^ 

Remark 3.5 Prom (3.6) one can easily see that if p and r] are both conserved densities of 
order n, with pn^n = P and rjn,n = Qi then ^ = hence the ratio of the top coefficients 
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is independent of x. If p and rj are conserved densities of consecutive orders say, \p\ = n and 
|r?| = n + 1 with pn,n = P and r?„+i,„+i = Q, then 



Remark 3.6 If the partial derivatives of F and p in (3.2) and (3.3) depend at most on Uj, then 
these equations are polynomial in Ujj^i, i > 0. In all the subsequent computations we have used 
only the coefficient of the top order derivatives. 



4. POLYNOMIALITY RESULTS FOR THE GENERAL CASE 



In this section we shall obtain polynomiality results, applying Proposition 3.3 either directly to 
a canonical density i = 1,2,3 (step 3 and 6), or to generic conserved densities p, v and 
show that one of the canonical densities is of that generic form (see Table 1). These derivations 
involve quite complicated and technical derivations organized in 6 steps. The intermediate results 
displaying polynomiality in u^, Um-i and Um-2 are presented respectively in Corollaries 4.1.2, 4.3.2 
and 4.6.2, labeled as "Result A", "Result B" and "Result C". 

Although the method seems to be a straightforward application of the equations (3.2) and 
(3.3) the computations are tedious and it took a number of trial and errors to discover the right 
sequence of computations presented as Steps 1 through 6. The computations are done in part 
analytically, in part with the symbolic programming language REDUCE. As the dependencies 
on lower order derivatives become explicit, the size of the actual expressions of the canonical 
densities grow very fast and we had to use generic forms for these conserved densities whenever 
possible. At each step, we tried to deduce from equations (3.2) and (3.3) a homogeneous linear 
system of equations, whose non-singularity leads to a polynomiality result. It is a remarkable 




hence Q = P- 
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fact in many cases these systems become singular for k = 2, hence our results are not applicable 
to fifth order equations. 

We use generic conserved densities in steps one, two and four. The first step is to obtain 
the quasi-linearity result for m > 5, which follows from the fact that the coefficient matrix of a 
homogeneous system is non-singular for m > 5. At the second and fourth steps we have a similar 
structure; we show that the coefficient of Um is independent of Um-i and Um-2 respectively, by 
obtaining nonsingular homogeneous systems of linear equations. 

The third and sixth steps are based on relatively straightforward computations using the 
canonical densities. At the third step we complete polynomiality in Um-i while at the fifth 
and sixth steps we complete polynomiality in Um-2 , by using the explicit form of the canonical 
densities p^^^ and p^^\ 

4.1. Step 1. quasi-linearity: = 

We now prove that an evolution equation admitting a quadratic conserved density is quasi- 
linear. The quasi-linearity result was already given in [5], but we shall repeat this derivation, 
because with Proposition 3.3, the proof given below is much shorter and neater. 

Proposition 4.1.1 Let ut = F{x,t,u, . . . ,Um), with m = 2A; -|- 1 > 13 and assume that it 
admits a conserved density of the form 

p = Pu^+i + QUm+l + R (4.1) 

where \P\ = \Q\ = \R\ = m. Then PFm,m = 0. 

Proof. Since the conserved density has order m+1, the restriction for the applicability of (A.6d) 
is A; -|- 1 > 7 hence m > 13. The coefficient of Um+i in (3.2) and the coefficient of Um+3 in (3.3) 
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are respectively as follows: 

i2k + 1) FmPm - {2k + 3) PFm,m = 

{2k + l){k^ + k + 6)FmPm - {2k + 3){k + l){k + 2)PFm,m = 
These two equations form the following homogeneous system of linear equations 



(4.2) 

(4.3) 



1 


-1 




F P 







{k'^ + k + 6) 


-{k"^ + 3k + 2) 




PF 

^ ^ m,m 








(4.4) 



The coefficient matrix in this linear system is singular for = 2. For k ^2 the homogeneous 
system has only the trivial solution hence from (4.4) we conclude that 



p p — pp 



0. 



(4.5) 



□ 

The quasi-linearity result follows from here immediately. 
Result A: Polynomiality in Um 

Corollary 4.1.2 If the canonical density p^^^ is a conserved quantity for the evolution equation 
lit = P, then Fm,m = 0) hence 

Ut = Aum + B 
where A, B are functions of x,t,u, . . . Um-\ ■ 

Proof, p^^^i in {A.2,) is of the form (4.1) with P = a'^a^ where a = fU"' ■ Hence PFm,m = 
implies Fm,m = 0; and ut = Aum + B □ 



In reference [5], by direct computation of the conserved density conditions, it was shown 
that equations of orders 7, 9 and 11 were also quasi- linear. Since the existence of the canonical 

13 



densities are necessary conditions for integrability, it follows that integrable evolution equations 
of order m = 2k+l with fc > 3 are quasi-linear. As discussed in Proposition 3.1, the applicability 
of the general formulas to the computation of the conserved density conditions require that 
k + I >7, hence the results below are valid only for m > 19. For lower orders we have checked 
the validity of our results by direct computations and outline the results in Section 5. 

4.2 Step 2: Polynomiality in Um-i, first result, Am-i = 

In the second and third steps we determine the dependency of the coefHcients A and B, in 
Uf = Aum + B on Ujn-i- For this purpose we consider a generic quadratic conserved density p 
of order m and use the coefficients of the top two nonlinearities in Dtp which are respectively 
given by equations (3.2) and (3.3). 

Proposition 4.2.1 Let ut = Aum + B, m > 19 with \A\ = \B\ = m — 1, and assume that it 
admits a conserved density 

p = Pul, + Qum + R, (4.6) 

where \P\ = \Q\ = \R\ = m — 1. Then 

PAm-i = 0. (4.7) 



Proof. We use here the equations (3.2) and (3.3) with I = 0, Fm = A and pn,n = 2P. The 
coefficient of Um in (3.2) is 



{2k + 1) Pm-iA - {2k + 3) Am-iP = 0, 



(4.8) 



while the coefficient of in (3-3) is 



2fe^ + 9A;2 + I3fe + 6 A^_iP - 2A;'' + 3A;^ + 13fe + 6 Pm-iA = Q 



3 I Q l2 



(4.9) 
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Prom (4.8) and (4.9) we get the following linear system. 



2k + 1 


-(2A; + 3) 











2P + 3fc2 + 13fc + 6 


-{2k^ + 9k^ + 13A; + 6) 




PAm-1 








(4.10) 



The coefficient matrix of the linear system is singular for A; = 2, then m = 5 which appear as 
an exception. For k ^2, the system has a trivial solution. Hence PA^-i = APm-i = 0. □ 
Again by using the canonical density p^^^ we shall obtain that A is independent of Um-i 

Corollary 4.2.2 // the canonical density p^^^ is conserved, then A^-i = 0. 

Proof. We substitute ut = Aum + -B in (^.3) and integrate by parts we can see that p*^^-* is of 
the form (4.6) where 

(4.11) 



a 



Hence PA^_i = implies Ajn_i = 0. 



□ 



4.3 Step 3: Polinomiality in m^-i, second result, -Sm-i,m-i,m-i — 

Now we shall see that the conservation of the canonical densities p^^^ and p^^^ will determine 
the form of B. By substituting ut = Aum + B, with |^| = m — 2 and |B| = m — 1 in the 
canonical densities p^^^ and p^^^ one can see that they reduce to 

where |QW| = = m - 2, and |Q(3)| = |i^(3)| = m - 1, and 

24 

P^^^ = -^^^—^a^-2,m-2 + 0^.20-1(^2 - 1), (4.12) 
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P(^) = la'io (m^ + Sm^ - 121m + 59?) 

m"^ + 3m^ — m — 3 L V / 

+ 60a-"'+^am-2Bm-i,m-i - l) 

+ 60a-2-+^i?^_,,_, (i^ - ^)] (4.13) 

Proposition 4.3.1 iei ut = Aum + B, m> 19 with \A\ = m — 2 and \B\ = m — 1. Then if 
the canonical densities p^^^ and p^^^ are conserved quantities, then 

Bm-l,m-l,m-l = 0. (4-14) 

Proof. To prove this result, it would be actually sufficient to use p^^^ only, but we first see how 
far one can go by using p^^\ We first compute (3.2) for p = p^^\ hence with Z = — 1, = A 

and Pn,n = Pm-\,m-\ = 2^^^^ to get 

{k + ^)APi^l^Um-l -{k- ^)P^^^Am-2Um-l - P^^^ Bm-1 = 0. (4.15) 

Differentiating (4.15) twice with respect to u^^i we obtain 

2P^^^Bm-l,m-l,m-l = 0. (4.16) 

If P(^) is nonzero, then B is quadratic in Ujn-i, but P(i) = gives a differential equation for 
a and we cannot exclude the possibility that Bjn-i,m-i,m-i / 0. It is possible to solve this 
differential equation, but it is easier to use p^^\ 

Now we compute (3.2) using p = p^^\ hence with I = 0, = A and pn,n = Pm,m = 2P(^) 
to get 

{k + ^)APi^l^Um -{k + i)p(3)^„_2u„_i - P(3)p„_i = 0. (4.17) 
This equation is linear in Um and its coefficient gives 

(fc + i)^Pil=0. (4.18) 
16 



Now since A should be nonzero, it follows that P^^^ is independent of Um-i- Differentiating 
p(^), in (4.13), with respect to Um-i, we obtain 



B. 



m—l,m—l,m—l 



am-2 - l) + 2a-2-+2B^_i,^_i (i- - -i,)] = 0. (4.19) 



Thus we can conclude that: 

-Bm-l,m-l,m-l = 0. (4.20) 

Result B: Polynomiality in Ujn,Ujn-i 

Corollary 4.3.2 // the canonical densities p^^^ and p^^^ are conserved quantities for the 
evolution equation ut = F, then 

F = Aum + C«^_i + Du„,^i + E. (4.21) 

where A, C, D, E are functions of x ,t, u, . . . Uni—2- 

4.4 Step 4: Polinomiality in w^_2, first result, A^_2 — C — 

At this step we assume the existence of two generic conserved densities p = Pu^_i + Qum-i+R 
and r] = Su^ + Tum + U and we use the relation S = Fm"^P given in Remark 3.5 and obtain 
CP = APfn-2 = APam-2/ci = 0. Then we compute the explicit form of the canonical densities 
p^^^ and p(^) to prove that these imply that Am-2 = C = 0. 



Proposition 4.4.1 Let ut = Aum + Cit^-i + Dum-i + E, m > 19, m = 2k + I with 
\A\ = \C\ = \D\ = \E\ = m — 2, and assume that it admits two conserved densities 

p = P«^_i + Qum-i + R, (4.22) 

ri = Sul, + Tum + U, (4.23) 
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where \P\ = \Q\ = \R\ = m - 2, \T\ = \U\ = m - 1 and S = F^'^P. Then 

CP = APm-2 = APam-2/a = 0. (4.24) 

Proof. If we compute (3.2) and (3.3) for Z = — 1 = A = a™, = ma"^~^am-2 = 

{2k + l)a"'-^ O'm—2 Pn,n — Pm—i,m—i — 2P we obtain the coefficient of Um—i in (3.2) as 

4CP - {2k + 1) APm-2 + {2k - l){2k + 1)AP^^ = (4.25) 

and the coefficient of Um+i in (3.3) as 

Uk'^CP - (2k^ + 3A;2 + 13A; + APm-2 

+ (2k^ - Zk^ + 13A; + e) {2k + 1)AP^^ = 0. (4.26) 

Then for Z = = A = and = pTn,m = 2S = 2a^P the coefficient of n^+i in (3.3) 
is obtained as 

I2{k^ + 2k + l)PC - {2k^ + 3k'^ + 25k + 12) APm-2 + 

+{4k^ + Ak^ + 23A;2 - 19A; - 6)AP^^ = (4.27) 

a 

Equations (4.25), (4.26) and (4.27) form the following system: 



4 


-(2A; + 1) 


(4P - 1) 




CP 









-{2k^ + 3k'^ + 13k + 6) 


{4k^ - 4A;3 + 23A;2 + 25A: + 6) 




APm-2 







12(P + 2A; + 1) 


-{2k^ + 3k^ + 25k + l2) 


{Ak^ + Ak^ + 23P - 19k - 6) 




1m-2 

a 








(4.28) 

The coefficient matrix in this linear system is singular for k = — ^, — |, 2. Notice that for k = 2, 
m = 5 is also an exception. For k ^ 2, the system has only the trivial solution. Thus from 
(4.28) we have (4.24) 
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□ 

Now we shall compute the canonical densities p^^\ p^'^\ and p^^^ and we shall see that their 
conservation will give ^^-2 = C = 0. 

By substituting ut = Aum + Cuf^_i + Dum-i + E, m > 19, m = 2k + 1 with A = a"^, 
\A\ = \C\ = \D\ = \E\ = m — 2, in the canonical densities p^^\ p^'^\ and p^^^ one can see that 
they reduce to 

p(2) = Sui,_^+Tul,_^ + Vum-i + U, 
with \Q\ = \R\ = \T\ = \V\ = \U\=m-2, \Y\ = \Z\ = m - 1, and 



P = a" + 



w? - 1 ia'f (1 - m) C , 
^ — + — a' 



24 



m 



1 C 2(m - 1) C2 



,2m 



(4.29) 



+ 
+ 
+ 



// / l(5 + m) C „ 
a a a a 

6 m a™ 

1 (m - 13) C, ^ (m - 

— a' a + 2^ ^ 

6 m Ou^ 

1 (6m - - 5) C 

6 

8 (1 
3 rrfi 



1) C 



a' a 



m 
m) 



2 C'C 



a 



3m 



(4.30) 



W = a{a'f- 12Q , , .7^.1. P7^— 

^ ' m{m^ + 3m2 - 121m + 597) oT- 

+ 240^^-^ -^o? (4.31) 
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where ' = _^ ■ Note that as higher order conserved densities should be quadratic (Corollary 
3.4) and top coefficients of conserved densities of consecutive orders are related (Remark 3.5), 
we have 

5 = and W = c?P. (4.32) 

We will now prove that, Proposition 4.4.1 together with the conservation requirements of the 
p(*^ 's above will imply A^n-i = C = 0. 

Corollary 4.4.2 Let ut = Aum + Cu^_i + Dum-i + E, m>19,m = 2k + l with A = a"^, 
\A\ = \C\ = \D\ = \E\ = m — 2. Then if p^^\ p^'^\ and p^^^ are conserved quantities, then 

Am-1 = C = 0. (4.33) 

Proof. Let's first consider the case where C = 0. In this case, P reduces to P = {a'^/a, and 
substituting this in (4.24), from APa'/a = we get a' = 0, hence Ajn-2 = 0. 

If C / 0, from (4.24), CP = implies P = 0. But as = a^P, we have W = also. Now, 
(4.29,4.30,4.31), is a system of three nonlinear differential equations for a and C. We can in 
principle solve a' from 1^ = in terms of C as the root of a quadratic equation, then substitute 
in P = and S = to get a system for C and C. With symbolic computations it is possible 
to see that C is zero. Wc give below an analytic proof of this fact. 

We define a = a' /a and C = C/a"^. Then C'/a"^ = C' + maC and a" /a = o' + a^. 
Hence 

P + l-2m~ . 
- = — — + aC + a' 

a 24: m 

+ ^-^^^C^-^C', (4.34) 



S o 1 (5m-m2 - 10)_2A 

D m 
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4m -2-2 ^im-lS)^~ 
D m 

6 m 

+ l^-^C', (4.35) 
3 m'* 



— = 0^-120 "^ + 3 

m{rin? + 3m^ — 121m + 597) 

^ ^^°m2(m3 + 3m2 - 121m + 597^ ' ^^'^'^^ 



Note that W/a^ is of the form 

W 



a 



^ = 0? - 2KiCa + K2C^ (4.37) 



where ki and K2 are constants that can be identified from (4.36) and since W = 



ai,2 = C{-Ki ± ^nl - H2). (4.38) 

If 7 denotes the coefficient of C in the expression above, we have a = and a' = 7C" and we 
obtain 



P / m2 - 23 o 1 - 2m m - 1 

- = 7' + 7 + 2 — 

a \ 24 m m^ 



+ [l--)c' = Q, (4.39) 
m/ 



~ 1 5 /o 7^ 4m + 5 4m -2 Sl-mA-o 

\ o m m^ 3 m'^ / 

+ (7' - -7 + A) = 0. (4.40) 
V m m^/ 

Now (4.39) and (4.40) form a system of homogeneous linear system for C' and C^. It can 
be checked that the coefficient matrix is nonsingular, hence (7 = 0. Prom (4.38), we have 
a = a' = 0, and the proof is complete. □ 
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4.5 Step 5. Polinomiality in second result: Djn-2,m-2 — 



We will now prove that D is linear in 'u^_2 
Proposition 4.5.1 Let 

ut = Aum + Dum-1 + E, m> 19, 



(4.41) 



with \A\ = m — 3, |D| 
form 



\E\ = m — 2, and assume that it admits a conserved density of the 



P = -P^m-l + QUm-1 + R, , 



(4.42) 



where \P\ = m — 3, \Q\ = \R\ = m — 2. Then 



D. 



m—2,m—2 



0. 



(4.43) 



Proof We compute (3.2) for equations (4.41) and (4.42) with / = -1, = A, pn,n = 2P to 
get: 



*4 



2Pm-3AUm-2 - 



" 2 



2PAm-3Um-2 = 2PD 



Differentiating (4.44) twice with respect to Um-2 we obtain 



PR 



m—2,m—2 



(4.44) 



(4.45) 



□ 



It can be checked that the canonical density p^^^ in (^.5) is of the form (4.42) and the 

coefficient P of the quadratic top derivative u^_i, is as follows: 



rn^- lOm^ + 720m - 2151 o 

P = ^ o — ««m-3 



+ 
+ 



- mm? + 9 
60(-3m2 + 8m + 3) 



m{m'^ — lOm^ + 9)a 
60 



m-3J^m-2 



m'^{m'^ + 4m + 3)a' 



2m" ^m-2 



(4.46) 
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That is 



P = pio) + pii) D' + p(^) D 



12 



(4.47) 



Hence PD" = implies 



(4.48) 



Prom this it follows that 



D" = n 



m— 2,m— 2 



(4.49) 



Corollary 4.5.2 // the canonical density p^^^ is conserved D is linear in Um-2- 



Thus we conclude that an integrable evolution equation should be of the form 



Ut = AUm + GUm-lUm-2 + Hu^-l + E, 



(4.50) 



where \A\ = \G\ = \H\ = m - 3 and \E\ = m - 2. 



4.6 Step 6. Polinomiality in Um-2, third result: 



'm— 2,m— 2,m— 2,m— 2 



The final step is to determine the Um-2 dependency in E using the explicit form of p^^^ . For 



Ut = AUm + GUm-lUm-2 + HUm-1 + E, m> 19, 



(4.51) 



with = |G| = \H\ = m — 3 and \E\ = m — 2, the canonical density p*-^^ is of the form 



■.(1) 



+ 



^m-2 

— 1 



2«m-2W,.m-3 



"1-3 /™2 



a.m-?,G , 1 



aG. 



K - 1) + 36^^(£- - 1) + 24:^ + 12-^(1 - £-) 



ma 



m 



ma 



2m 



m 



m^ — 1 
+ 12u^_2 



(m -1) + 18— -— (— -1) + 12 



a'" m 



ma" 



a'^"^m'^{m + 1) 



2 '^m-4 



m-3" 



12u 



[2aG - ma^am-3] 



a 



m— 3 



+ 12 



m^im? — 1)0^'^ 



m{m + l)a"^ 

2ma'^Em-2 + H'^{m-l) 



(4.52) 
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where a"* = A 
Proposition 4.6.1 Let 

Ut = AUm + GUm-lUm-2 + HUm-1 + E, m> 19, (4.53) 
with \A\ = \G\ = \H\ = m — 3 and \E\ = m — 2, Then if p^^^ is a conserved quantity, then 

Em-2,m-2,m-2,m-2 = 0. (4.54) 

Proof. Notice that p^^^ has order m — 2 but it is not a priori even polynomial in Um-2 since 
we don't know the form of E. We substitute (4.53) in (3.2) with I = —2, = A and 
Fm-i = Gum-2 + H 

{k + l)ADp^^l,^^_^ -{k- l)DApil^^^_^ = {GUm-2 + i^)pS-2,n^-2 (4-55) 
The coefHcient of in (4.55) is: 

P^Um-2,m-2 = 0, (4-56) 



hence since a ^ 0, it follows that 



dum-2^ m{m? — 1) 



Em-2,m-2,m-2,m-2 = 0. (4-57) 



□ 



Result C: Polynomiality in Um,Um-i,Um-2 

Corollary 4.6.2 // the canonical densities p^^\ p^'^^ and p^^^ are conserved quantities for the 
evolution equation ut = F , then 

F = AUm + GUm-2Um-l + HUm-l + Jul^-2 + L'>4a-2 + NUm-2 + S (4.58) 

where A, G, H, J, L, N, S are functions of x,t,u, . . . 11^^—3. 
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Remark 4.7 By the use of the conserved densities p^'^\ z = 1,2, 3 we have shown that A is 
independent of iim-3 and G = 0. However it was not possible to obtain polynomiality in 
because further conserved densities were needed. 

5. CONCLUSION 

In the present paper we obtained ah the polynomiahty information that could be extracted 
from the conserved densities up to p^^^ (See Remark 4.7). In general computations we haven't 
used p(~^) and p^^\ but using these for m = 7 didn't give any further information. 

For m = 7, in initial computations we used dependencies in all derivatives. But at later 
stages this was impossible and in the search of computationally efficient methods, we noticed 
that all polynomiality results followed from the coefficient of the top order term in (3.2) and 
(3.3). The use of only the top dependency lead us to formulate a graded algebra structure on 
the polynomials in the derivatives Uk+i, that we called "level grading" [8]. This structure is 
based on the fact that derivatives of a function oi x,t,u, . . . , Uk are polynomial in higher order 
derivatives and have a natural scaling by the order of differentiation above the "base level k." 
The crucial point is that equations relevant for obtaining polynomiality results involve only 
the terms with top scaling weight with respect to level grading. Thus one can work with the 
dependency on the top level term only and reduce the scope of symbolic computations to a 
feasible range. Applications of the "level grading" structure to the classification problem will 
be discussed elsewhere. 

By the remark following Proposition 3.1, the polynomiality results obtained in Section 4 are 
valid for m > 19. For m = 7, 9, 11, 13, 15, 17 the conserved density conditions are computed, 
directly, without using (3.2) and (3.3), with the symbolic programming language REDUCE and 
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it is shown that the evolution equation have again the form given by (4.58). The results of 
Corollary 4.6.2 are valid for all m > 7 and we restate it here for convenience. 

Corollary 4.6.2 Let 

Ut = F{x,t,U,... ,Um) 

be a scalar evolution equation in one space dimension of order m = 2k + l where m> 7. If the 
canonical densities p^^\ i = 1,2,3 given in the Appendix are conserved quantities then 

F = AUm + GUm-2Um-l + HUm-1 + JUm-2 + + + S. 

where A, G, H , J , L, N , S are functions of x,t,u, . . . Um-3- 

For future work, the graded algebra structure seems to be promising tool both for the 
computation of conserved densities and for obtaining integrability conditions. In fact, conserved 
densities for seventh order equations have been computed and it has been observed that there 
are candidates for integrable equations where A may depend on U4 and is non-polynomial in U4. 
In this sense. Corollary 4.6.2 is the best polynomiality result for evolution equations of order 7. 



APPENDIX A: Canonical Densities of the Formal Symmetry Method 

If the evolution equation ut = F[u] is integrable, it is known that the quantities 

p(-i)=F-V-, p(o)=F^_i/F^, (Ala) 

where 

dF dF 
Fm = ^, F^_i = - (A.lb) 

dUm OUm-1 

are conserved densities for equations of any order [3]. The next three canonical densities up to 
total derivatives computed in [5] are presented below. The expressions of the canonical densities 
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up to total derivatives are convenient for the present work but one has to be careful in general 
in adding a total derivative to the canonical densities. Because if p^*) is a canonical density with 
DtP^"^^ = Da^^\ then p(*+'=) will involve cr^*) after some k [3]. Hence although total derivatives 
in are irrelevant for the existence of cr^*^ , they should not be omitted whenever cr^*-' enters 
in the expression of another canonical density. We shall use the following notation 



F 



, z = 1,2,3,4. 



(^.2) 



= a-\Daf- 



12 



m(m + 1) 



Daa 



(i)+a 



12 2 24 
m?{m + 1) m,{m? — 1) 



a 



(2) 



(A.3) 



a{Da) 



Dam H "(n 



+ 2a^ 



2«(1) + 



m(m — 1) 



(m- 1) 



a(i)«(2) 



a 



(2) 



(m — 1)(to — 2) 



a 



(3) 



(^.4) 



9(3) = aiD-'af- 

+ 30a{Daf 
120 



+ 



+ 



60 1 
m(m+l)(m + 3) ^^^4 ^ ' 

(m- 1) 



/ i\/ n\"-^'^(l) "I o7 7 — o T\"'-''f2') 

m(m + l)(m + 3) ^ m^(m + l) ^ m(m^ — 1) 



m{m^ — l)(m + 3) 



(m — l)(m — 3) 



a(x)Da(x) + (m — 3)Da 



(2) 



(m-l)(2m-3) 3 6(m - 2) 



60 



«(!) + 



m{m? — l)(m + 3) 



m 

{m — 1) 



-«(i)a(2) - 6q;(3) 
4 



/T-v \2 - r-v (m — l)(2m — 3) 4 
(^a(i)) Daa)a(2) H 3 



- 4 



(2m -3) 2 



-7;^«(i)«(2) + -«(i)«(3) + -«(2) - 



(4) 



(A.5) 
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Evolution equation of order m with generic and canonical densities 



Step 


Evolution Equation 


Type of the 
conserved density 


Form of the 
conserved density 


Result 


1 


Ut = F[u] 


generic 


P = Pu^+i + QUm+l 


+ R 


OF _ 


-- 


2 


Ut = Aum + B 


generic 


p = Pul^ + QUm + R 




dA 


= 




dUm-l 


3 


Ut = Aum + B 


canonical 


p = p(^) 






= 






4 


Ut = AUm + Cu^_i 
+DUm-l + E 


generic and 
canonical 


r] = Sul, + Tum + U 


+ R 


dA 


= 0, 


dUm-2 

c = o 


5 


Ut = AUm + DUm-1 

+E 


generic and 
canonical 


p = Pu^rn-l + QUm-l 


+ R 




= 




6 


Ut = AUm 

+GUm-lUm-2 
+HUm-l + E 


canonical 


P = P^'^ 






= 




a<-2 



Tablel:Polynomiality results for the general case 
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APPENDIX B: Expressions of fc'th Order Derivatives [5] 



D^ip = (pnUn+k + P{Un+k-l), k>l 



{A.6a) 



{AM) 



n\ ^n+fc— 1 



+ 



iPn-2 + kDipn-1 + I -D (/?n 



Un+k-2 + PiUn+k-s), k>5 



{A.6c) 



D (p = (fnUn+k + i'Pn-l + kDtfn] Un+k-l 



+ 



+ 



(pn-2 + kDipn-1 + {^\D (pn 



Un+k-2 



Un+k-3 



+ P{Un+k-4:), k>7 



(AM) 
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